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ABSTRACT 

We present a fully consistent evolutionary disc model of the solar cylinder. The model is 
based on a sequence of stellar sub-populations described by the star formation history 
(SFR) and the dynamical heating law (given by the age-velocity dispersion relation 
AVR) . The stellar sub-populations are in dynamical equilibrium and the gravitational 
potential is calculated self-consistently including the influence of the dark matter halo 
and the gas component. The combination of kinematic data from Hipparcos and the 
finite lifetimes of main sequence (MS) stars enables us to determine the detailed ver- 
tical disc structure independent of individual stellar ages and only weakly dependent 
on the IMF. The disc parameters are determined by applying a sophisticated best fit 
algorithm to the MS star velocity distribution functions in magnitude bins. We find 
that the AVR is well constrained by the local kinematics, whereas for the SFR the 
allowed range is larger. The model is consistent with the local kinematics of main 
sequence stars and fulfils the known constraints on scale heights, surface densities and 
mass ratios. A simple chemical enrichment model is included in order to fit the local 
metallicity distribution of G dwarfs. In our favoured model A the power law index of 
the AVR is 0.375 with a minimum and maximum velocity dispersion of 5.1km/s and 
25.0 km/s, respectively. The SFR shows a maximum 10 Gyr ago and declines by a 
factor of four to the present day value of 1.5MQ/pc^/Gyr. A best fit of the IMF leads 
to power-law indices of -1.46 below and -4.16 above 1.72Mq avoiding a kink at IMq. 
An isothermal thick disc component with local density of ~ 6% of the stellar density 
is included. A thick disc containing more than 10% of local stellar mass is inconsistent 
with the local kinematics of K and M dwarfs. Neglecting the thick disc component 
results in slight variations of the thin disc properties, but has a negligible influence 
on the AVR and the normalised SFR. The model allows detailed predictions of the 
density, age, metallicity and velocity distribution functions of MS stars as a function 
of height above the mid-plane. The complexity of the model does not allow to rule 
out other star formation scenarios using the local data alone. The incorporation of 
multi-band star count and kinematic data of larger samples in the near future will 
improve the determination of the disc structure and evolution significantly. 

Key words: Galaxy: solar neighbourhood - Galaxy: disk - Galaxy: structure - 
Galaxy: evolution ~ Galaxy: stellar content - Galaxy: kinematics and dynamics 



1 INTRODUCTION 

There is no 'concordance' Milky Way model available so far 
that describes the structure, kinematics, chemistry and evo- 
lution of the disc in great detail . The c las sical stellar den- 
sity model of lBahcall fc Soneiral (Il980al lbl Il984h composed 
by a spheroid and an exponential disc with magnitude- 
dependent exp onential scale heights is still widely used. 
iBahcalj (|l984al lbh introduced a finite set of isothermal disc 



components and solved the Poisson and the Jeans equation 
for a dynamical equilibrium model of the vertical disc struc- 
ture. At present the most sophisticated Milky Way model 
based on star counts is the so-called 'Besangon model' of 
th e Galaxy developed and presented in a series of articles. 
In lRobin et al.l l|2003l ) a general description and the present 
status of the model is given. In this model the thin disc is 
composed by a sequence of stellar sub-populations with in- 
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creasing age and velocity dispersion. The Besangon model 
has still some serious drawbacks in the construction of the 
thin disc concerning the density profiles of the components, 
the star formation history (SFR) and the initial mass func- 
tion (IMF) that will be discussed below. 

The main input functions to be determined in an evo- 
lutionary disc model are the SFR and the dynamical heat- 
ing described by the age- velocity dispersion relation ( AVR) . 
The SFR(t) of the Milky Way disc is still not very weU 
determined. The main reason for that is the lack of good 
age estimators with corresponding unbiased stellar samples. 
Especially samples selected by colour cuts or by a mag- 
nitude limit are biased with respect to the age distribu- 
tion, because there is an age-metallicity relation due to 
the chemical enrichment of the disc. Therefore not even 
the famous Geneva-Copenhagen sample of 14,000 F and G 
stars iNordstrom ct al. 2004. ; Holmberg et al. 2009) (here- 
after GCS1,GCS2) with well determined stellar properties 
and individual age estimates can be used to derive the star 
formation history directly by star counts. Systematic biases 
introduced by the method of GCSl concerning age and stel- 
lar pa ramet er determinations are discussed in Po nt fc Eveil 
l|2005h and iHavwoodl (|2006l '). The Besangon model of the 
Galaxy was dev e loped a nd presented in a series of articles. In 
iHavwood et al.1 (|l997al lbh the SFR of the thin disc is deter- 
mined to be approximately constant using a series of isother- 
mal components and applying an approximate method to 
achieve dynamical equilibrium. In further applications they 
used a constant SFR combin ed with a steep IMF in the sensi- 
tive mass regime of 1-3 Mq (jRobin et al. 2003). As a result, 
the mean age of the disc population, the scale heights and 
the surface d ensities of the stellar sub -populations are rel- 
atively small. iRocca- Pinto et al.l (|200G| ) used chromospheric 
age determinations of late- type dwarfs. They apply stellar 
evolution and scale height corrections. The main result is 
the determination of enhanced star formation episodes over 
the lifetime of the disc. They exclude a constant SFR from 
chemi cal evolution models. In lR. fc C. de la Fuente Marcos! 
l|2004h the star formation history for open star clusters was 
determined. They found at least five episodes of enhanced 
star formation in the last 2 Gyr. Since star clusters contain 
only a small percentage of all disc stars, an extrapolation to 
the total (smoothed) SFR is not possible. 

In recent years the Hipparcos stars with precise paral- 
laxes and proper motio ns were used to determ ine the SFR 
with different methods. iHernandez et all (|2000l ) determined 
the SFR over the last 3 Gyr using isochrone ages. They 
found a series of star formation episodes on top of an un- 
derlying smooth SFR. This result is complementary to our 
model, which g ives t he slow evolution of the smoothed SFR. 
IVereelv et aP \2m± applied an inverse method to derive 
from the colour magnitude diagram the star formation his- 
tory and age-metallicity relation. They used prescribed AVR 
and gravitational potential and fixed the IMF to be a power 
law over the whole stellar mass range. This ansatz results 
in a steep IMF and a rather young thin disc. Addition- 
ally they found s trong maxima in the SFR at ag es of 10^ 
and 2 x lO'^ yr. In lBinnev et al.l ([2000,) and in .Cienoni et al.l 
l|2006t ) the scale height variation of main sequence stars was 
not taken into account. Therefore these models derive the 
local age distribution of stellar sub-populations with life- 
times larger than the age of the disc instead of the SFR. 



iBinnev et all (|2000l ) determined the age of the solar neigh- 
bourhood to be 11.2 ± 0.75 Gyr, which includes the con- 
tribution of the thick disc. Both papers found an approx- 
imately constant local age distribution, which corresponds 
to a decreasing SFR according to the increasing scale height 
with incre asing age of the stellar su b-populations. In a re- 
cent paper lAumer fc Binnevi ||2009|') derived the SFR using 
the updated Hipparcos data (|van Leeuwen ^ 2007) and the 
kinematics of the GCS sample. They calculated in detail 
the number of main sequence (MS) stars as function of B-V 
colour based on a KTG93-like IMF and scale height correc- 
tions due to dynamical heating. Their favoured disc model 
has an exponentially decaying SFR with a decay timescale 
of 8.5 Gyr and an age of 12.5 Gyr. In the model the gravi- 
tational potential is fixed and the SFR depends strongly on 
the properties of the IMF by construction. 

The presented paper (Paper I) starts a new attempt to 
develop a fully consistent disc model which is able to predict 
the properties of the Milky Way disc in great detail concern- 
ing density distributions, age distributions, kinematics and 
chemistry of all types of stars. The aim of this sequence of 
papers is to construct a smooth physically consistent disc 
model that can be extended to a global Galaxy model and 
allows detailed predictions of number densities and kine- 
matics of stars of different types. The new model can be 
of some help for the construction of a 'concordance' model 
of the Milky Way. It can be used for the preparation of 
the future astrometric space mission Gaia that will measure 
with high precision positions, distances, proper motions and 
stellar parameters (temperature, surf ace gravity and chem- 
ical co mpositio n) o f one billion stars (jPerrvman et al.ll200ll : 
iBailer- Jones .20051 ). 

We start with a local disc model of the solar cylin- 
der, which is based on the SFR, the AVR and a chemical 
evolution function described by the age-metallicity-relation 
(AMR). The vertical density profiles and the correspond- 
ing scale height determinations of thin and thick disc stars 
are calculated in a self-consistent gravitational potential in- 
cluding the DM halo and the gas component. The model 
parameters are determined by a best fit procedure of the 
velocity distribution functions /i(|W|) of solar neighbour- 
hood stars from B-K dwarfs. We use Hipparcos stars and at 
the faint end the Catalogue of Nearby Stars (CNS4), select 
the MS stars and divide these into a series of volume com- 
plete sub-samples in magnitude bins. The derived /i(|V[^|) of 
each sub-sample are compared simultaneously to the model 
predictions. Details of the fitting procedure and the signifi- 
cance with resp ect to parameter var iations are discussed in 
a future paper IJust fc Jahreidl2009l ) (Paper II). 

Since the velocity distribution function of MS samples is 
the average over the lifetime of these stars properly weighted 
by the SFR and the dilution due to the increasing scale 
height, the time resolution is strongly limited. Therefore we 
use only smooth input functions for the SFR, AVR and 
metal enrichment in the model. The resulting SFR, AVR 
and chemical evolution describe the long-term smoothed disc 
evolution. Similar disc models by fitting the vertical lumi- 
nosity and colour profiles of edge-on galaxies ins tead of the 
kinematics have already been used successfully ijjust et al.l 
Il996l . l2006l ). 

The main advantages of this method is that it does not 
depend on individual stellar age determinations and that 
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it is essentially independent of the shape of the IMF. A 
weakness is that we cannot confine the SFR strongly based 
on local data only. The application of the model to star 
counts of remote stellar populations will be of great help 
in this respect. In subsequent work the local model will 
be continuously extended and compared to large data sam- 
ples taken fi-om the Slo an Digital Sky Survey SDSS/SEGUE 
llAbazaiian et al] 2009) and the Radial Velocity Experiment 
RAVE (jZwitter et al. 2008,') to further constrain the possi- 
ble parameter range of the disc model. 

A major restriction of a local model is that radial mix- 
ing of stars in the disc cannot be included easily. There 
are two main mixing processes discussed in the literature. 
The first one is directly connected to the dynamical heating 
process responsible for the AVR. The gravitational scatter- 
ing process leads to a diffusion of orbits in velocity space 
(AVR) and in posit ion (radial, tangential, vertical mixing). 
IWielen et al.l (j 19961 ) discussed the radial diffusion of stellar 
orbits and the consequences for the AMR in some detail. The 
radial probability distribution function of the birth-places 
depends mainly on the age of the stellar sub-population and 
only weakly on the physical scattering process. For a 5 Gyr 
old population the typical radial width of the initial distribu- 
tion is ±2 kpc. The second mecha nism is resonant scatterini 
of circular orbits by spiral arms llSellwood fc Binnevl[200_. 
iRoskar et all l2008l : ISchonrich fc BinnevI l2009t ) . This mech 
anism changes the radial position of stars quickly, but the 
eccentricity of the orbits remain very small. The effect of res- 
onant scattering may account for up to h alf the stars in the 
solar neighborhood (jRoskar et al.ll2008h . which shows the 
possible scale of the errors which may result if migration is 
ignored. Unfortunately the efficiency and the properties of 
resonant scattering in the Milky Way disc are still poorly 
known. In a future global disc model radial mixing may be 
included in a parametrized form. The local model presented 
in this paper is primarily a description of the 'status quo' 
of the local Milky Way disc quantified by the local age dis- 
tribution (SFR), the kinematics (AVR) and the chemical 
composition (AMR). Therefore the relations between these 
function are not altered by radial mixing. But the physi- 
cal interpretation of the SFR as a 'local star formation his- 
tory' and the AVR as a 'local dynamical heating process' 
are weakened by radial mixing and require corrections. For 
the chemical evolution model the consequences are more se- 
vere. The tight connection of the enrichment, the SFR and 
the gas infall will be broken. The chemical enrichment of 
the gas decouples partly from the (apparent) enrichment of 
the stellar population. A more general local AMR including 
an intrinsic scatter would b e the consequence as derived by 
ISchonrich fc BinnevI (|2009l ). 

In section [5] the physics of the disc model is presented, 
section [3] describes the observational data, in section |3] the 
fitting procedure and the properties of the best fit models 
are given, section [S] collects the main results and discussion. 



2 DISC MODEL 

In this section we describe the construction of the disc model 
and the determination of the properties of the stellar sub- 
samples. We use a thin, self-gravitating disc composed of a 
sequence of stellar sub-populations according to the SFR 



and the AVR as input functions. Additionally the gravi- 
tational force of the gas component and the DM halo are 
included. The bulge and the stellar halo can be neglected 
in the solar neighbourhood. The total gravitational poten- 
tial and the density profile of each age-bin are determined 
self-consistently assuming isothermal distribution functions 
with velocity dispersion according to the AVR. We take into 
account finite stellar lifetimes and mass loss due to stellar 
evolution. The stellar remnants stay in their sub-population, 
the expelled gas (stellar winds, PNs, and SNs) is mixed im- 
plicitly into the gas component. Since the lifetimes and mass 
loss depend on metallicity, the metal enrichment with time 
is included. A standard IMF is used for the determination 
of the stellar mass fraction of the sub-populations with age. 
The velocity distribution functions /ms(|M^!) for main se- 
quence stars are calculated by a superposition of the Gaus- 
sians weighted by the local density contribution up to the 
lifetime of the stars. 



2.1 Self-gravitating disc 

The backbone of the disc is a self-gravitating vertical disc 
profile in the thin disc approximation including the gas com- 
ponent and the DM halo. In this approximation the Poisson- 
Equation is one-dimensional and in the case of a purely self- 
gravitating thin disc (i.e. no external potential) the Poisson 
equation can be integrated leading to 



d£ 

dl 



SttG / p($')d$' 



(1) 



Therefore we model all gravitational components by a thin 
disc approximation. We include in the total potential $ the 
contributions of the stellar component <I>s, the gas compo- 
nent <I>g, and the DM halo (E>h 



^{z) = $4^) + $3(2) + $1,(2) 



(2) 



In order to obtain the force of a spherical halo correctly in 
the thin disc approximation, we use a special approximation 
(see subsection 12. 4|l . Since we construct the disc in dynami- 
cal equilibrium, the density of the sub-components are given 
as a function of the total potential Pj(4'). The vertical dis- 
tribution is given by the implicit function z('&) via direct 
integration 



f d$' SttG / 
Jo Jo 



p($")d$" 



-1/2 



(3) 



For the fitting procedure it is essential to separate the nor- 
malised model and the scaling factors. Since we use the SFR 
and the AVR as input functions, it is convenient to nor- 
malise the model to the total surface density Etot and to 
the maximum velocity dispersion cto of the oldest thin disc 
sub-population. As a consequence the gravitational poten- 
tial in equations [T] and [3] is normalised to and a natural 
scale length is 



(4) 



27rGEt, 



which corresponds to the exponential scale height of an 
isothermal component with velocity dispersion (Je above a 
gravitating sheet with total surface density Etot. 
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The normalised disc model describes the intrinsic struc- 
ture of the disc. The relative thin disc, thick disc, gaseous, 
and DM fractions of the surface density (up to \z\ = Zmax) 
are given by the input parameters Qs,Qt, Qg, Qh that have 
to be iterated to find the best model. The value of Zmax is 
determined by the prescribed maximum of the normalised 
potential $(«max)/o"c- The disc model has two free global 
scaling parameters to convert the normalised model to phys- 
ical quantities. One can choose two parameters out of four 
types: a local density, a scale height or thickness, a surface 
density, or a velocity dispersion (see section |4ll. 



2.2 Stellar disc 

The stellar component is composed of a sequence of isother- 
mal sub-populations characterised by the IMF, the chemi- 
cal enrichment [Fe/H](t), the star formation history SFR(t), 
and the dynamical evolution described by the vertical veloc- 
ity dispersion (tw(t). Throughout the paper we use t for the 
time with present time tp — 12 Gyr and r = tp — t for the 
age of the sub-populations. We include mass loss due to stel- 
lar evolution and retain the stellar-dynamical mass fraction 
(?(r) (stars + remnants) only. The mass lost by stellar winds, 
supernovae and planetary nebulae is mixed implicitly to the 
gas component. 

With the Jeans equation the vertical distribution of 
each isothermal sub-population is given by 



Psj(^) = Psojexp I — 2 



(5) 



where ps,j is actually a 'density rate', the density per age bin, 
and o"w,j the velocity dispersion at age rj. The connection 
to the SFR is given by the integral over z 



(6) 



with time tj = tp — Tj. The (half-)thickness /lagc(Tj) of the 
sub-populations is defined by the mid-plane density psoj 
through 



ff(rj)SFR(tj) 

2/lago(''j) 



and can be calculated by 



/lagc(Tj) — 



pB.j(z) 



dz = 



exp 



-<S>{z) 



'w,j 



dz 



(7) 



(8) 



The total stellar density Ps{z) and velocity dispersion as{z) 
are determined by 



O-s(z) = 



ps.i{z)dt 



P^{z) Jo 



'^w{Ti)ps,j{z)dt 



(9) 

(10) 



The stellar surface density Eg, which includes luminous stars 
and stellar remnants, is related to the central density ps,o 
by the (half-)thickness h^s and can also be converted to 
the integrated star formation So using the effective stellar- 



dynamical mass fraction QcB 



So 



PB{z)dz = 2/icffPsO = ffcffSo 



SFR(t)dt = {SFR)tp 



(11) 
(12) 

(13) 



ffcft = ^ / giT)SFR{t)dt 
So Jo 

The general shape of pa{z) can be characterised by two 
shape parameters which are h^B and the exponential scale 
height Zs well above the plane. For an exponential profile 
fteff = Zs and for an isolated isothermal sheet oc sech^ {z/2zb) 
we find hcH — 2zs. Most realistic profiles are somewhere 
in-between. From the normalized version of equation [TT] we 
calculate hcs/zc- The effective exponential scale height Zs/z^ 
of the stellar disc is determined numerically by the mean 
exponential scale height in the range z/zc = 2 ... 5. 

The metallicity [Fe/H] affects the lifetimes, luminosities 
and colours of the stars and as a consequence also the mass 
loss of the sub-populations. In order to account for the influ- 
ence of the metal enrichment we include a metal enrichment 
law AMR that leads to a local metallicity distribution of late 
G dwarfs consistent with the observations. The properties of 
the stars and the stellar sub-populations are determined by 
population synthesis calculations (see Sect. 12.5)1 . 

The properties of MS stars with lifetime r are deter- 
mined by an appropriate weighted average over the age 
range. We use 



Pr(z) 



SFR(t')dr' 



SFR(t') 

2/lagc ("^ 



exp 



dr' 



(14) 
(15) 



where Sr are all stars born up to age t and Pt{z) would 
be the density profile of these stars, if mass loss by stellar 
evolution is ignored. The thickness hms and the normalised 
density profile pmsiz) / Pms{0) of MS stars are given by 



2pr{0) 



and 



Pms(z) _ Pt{z) 
Pms(O) ^ p,(0) 



(16) 



The normalised density profile is independent of the IMF. 
For the determination of the number density profile of MS 
stars with lifetime t the fraction of stars AA'^ in the cor- 
responding mass interval of the IMF must be taken into 
account (similar for the mass density profile). 

In a similar way we calculate the velocity dispersion pro- 
file ains{z) and the velocity distribution functions fnis{W, z) 
for the MS stars using 



U.{w,z) = 



Pr{z) 
1 



'^w(T')PT'(z)dr' 



Pt'[z) 



Pr[z) Ja aw{T')V2^ 



exp 



2<(r') 



(17) 

dr' (18) 



For the thick disc component we use a simple isother- 
mal component adopting an age larger than tp. The thick 
disc is parametrised by the surface density Et and the ve- 
locity dispersion at. The density profile is calculated self- 
consistently as for the thin disc sub-populations. The thick 
disc contributes only to the velocity distribution functions 
of the lower MS with lifetimes larger than tp. For the rel- 
ative contribution to the fi{W) we add a correction factor 
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(~ 1.8) to the local mass ratio ol thick and thin disc to ac- 
count for the different age distribution of thin and thick disc. 
The larger mean age of the thick disc leads to an enhanced 
number density of low mass stars due to the larger mass loss 
by stellar evolution. 



The DM halo is parametrised as an isothermal compo- 
nent similar to the thick disc by the velocity dispersion a"h 
and the surface density Eh, which is implicitly determined 
by the total surface density. 



2.3 Gas component 

For the gas component we use a simple model to account 
for the gravitational potential of the gas. The vertical pro- 
file of the gas component that is used for the gravitational 
force of the gas, is constructed dynamically like the stellar 
component. The gas distribution is modelled by distribut- 
ing the gas with a constant rate up to a maximum age over 
the velocity dispersion range Sgcrw(r) of the young stars 
reduced by a factor Sg to account for a smaller minimum 
value. By varying these parameters the peakyness and the 
width of the gas density profile can be changed. We adjust 
the (half-)thickness and the exponential scale-height of the 
gas compon ent to the observed values of /loff.g ~ 150 pc and 
Zg ~ 100 pc iDickev fc Lockmanlll99(il ). We determine % nu- 
merically in a similar way as Zs but closer to the mid-plane 
in the range of ^/(O.S z^) = 2 ... 5. The surface density Eg 
of the gas is iterated to reproduce the observed value. 



2.4 Dark matter halo 

The halo does not fulfil the thin disc approximation. For a 
spherical halo we get the vertical component of the force to 
lowest order from 



d$h 
dz 



i?2 



(19) 



with — E?+z^ and circular velocity iic,h of the DM halo at 
radius R. Comparing this with the one-dimensional Poisson 
equation from the thin disc approximation (integrated over 
z near the mid-plane to lowest order for small z) 



d£ 
dl 



AirGpoz 



we can use for the local halo density 



«c.h 

47rGi?2 



(20) 



(21) 



to be consistent with the spherical distribution. Only for 
a singular isothermal sphere this virtual halo density cor- 
responds exactly to the real local DM density. The second 
parameter is the halo velocity dispersion (Jh . Comparing the 
spherical distribution of the halo density and the vertical 
profile of a self-gravitating isothermal sheet up to second 
order in z we find the relation 



dlnph/dlni? 



(22) 



For an isothermal sphere the logarithmic derivative of the 
density equals two. As a consequence we can apply the thin 
disc approximation also for the halo by using the virtual 
halo density pho and velocity dispersion (Th estimated from 
the circular velocity Vc.h- For a singular isothermal sphere 
these quantities correspond to the physical values. The effect 
of different radial profiles, anisotropy and flattening would 
lead to correction factors. 



2.5 Stellar population synthesis 

For the determination of luminosities, main sequence life- 
times and mass loss we are not interpolating directly evo- 
lutionary tracks of a set of stellar masses and metallicities. 
In order to get a complete coverage of stellar masses we 
u se instead the stellar populatio n synthesis code PEGASE 
(Fioc fc Rocca-VolmeranggI 19971 ) to calculate mass loss and 
luminosities of "pseudo" simple stellar populations (SSP). 
This means that the PEGASE code is used to calculate the 
integrated luminosities and colour indices for a stellar popu- 
lation created in a single star-burst as function of age. These 
SSPs are then used to assemble a stellar population with a 
given star formation history, in the sense that the star forma- 
tion history is assembled by a series of star-bursts. In this 
way we can assemble stellar populations for varying SFR 
and metal enrichment [Fe/H]. Our SSPs are modelled by 
a constant star formation rate with a duration of 25 Myr, 
the time resolution of the disc model. This is done for a set 
of different metallicities and intermediate values from the 
chemical enrichment are modelled by linear interpolation. 

The application of the PEGASE code is twofold. On the 
one hand mass loss due to stellar winds, planetary nebulae 
and supernovae determines (/(r), the mass fraction remain- 
ing in the stellar component as a function of age r. This 
depends o n the IMF and the metallicity. We adopt a Scalo- 
like IMF (|Scalolll986l ) given by 



diV 



M°dM 




0.08 M/Mq < 1 
for 1 M/Mq < 2 
2 s; M/Mq < 100 



(23) 
(24) 



and five different metallicities ([Fe/H]= -1.23, -0.68, -0.38, 
0.0, 0.32), which are the input parameters of the code. Fig.[T] 
shows the mass loss for the different metallicities. The frac- 
tions of luminous matter and of stellar remnants for the set 
of input metallicities are shown separately. The sum of both 
contributions for each metallicity vary by a few percent only 
and are therefore not shown. The thick full (red) line shows 
for the fiducial model A the fraction of stellar mass g{T) in 
the present day stellar disc as a function of age including the 
chemical evolution. For the oldest age-bins about 40% of the 
stellar mass is lost by stellar evolution. The total fraction of 
stellar mass g^s = 0.654 is indicated by the horizontal line. 

In the second application we determine the MS lifetimes 
and luminosities for stellar mass bins. Here we use the PE- 
GASE code to compute V-band luminosities of isochrones 
in small mass bins. These are needed to estimate the mean 
MS lifetime as a function of the absolute V-band luminos- 
ity Afv for the calculation of the velocity distribution func- 
tions f{\W\) of these stars. We apply the PEGASE code 
to piecewise constant IMFs for mass- bins with O.IMq and 
for a finer grid of metallicities ([Fe/H]= -0.8, -0.68, -0.5, - 
0.28, 0.0, 0.20, 0.32). In Fig.[2]we show Mv for the different 
mass bins as a function of age. The lower panel gives the 
luminosity evolution of a M = 0.8 Mq star for the differ- 
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Figure 1. Mass loss due to stellar evolution. The two sets of 
thin lines show the fractions of luminous mass and of remnants 
(upper and lower curves, respectively) for the set of metallicities 
used in the PEGASE code. The full (red) thick line shows the 
total stellar mass fraction (luminous matter plus remnants) of the 
fiducial model A taking into account the age-metallicity relation. 
The overall present day mass fraction of the stars is given by g^g . 



ent metallicities (thin lines) and the age dependence of the 
luminosity for different disc models for M — 0.8 Mq demon- 
strating the significant variation over the whole age range. 
The upper panels show the age dependent luminosities for 
all mass bins. These are again not the luminosity evolutions 
of the stars but the present day properties of the stellar pop- 
ulation taking into account the age-metallicity relation. The 
vertical lines indicate the estimated mean MS lifetimes in 
the corresponding luminosity bins. 

For a consistent model the locations of the mass tracks 
(B — V,Mv)(t) in the Hertzsprung- Russell diagram (HR- 
diagram) should cover the same area used for the selection of 
main sequence stars (see Fig. |3} ■ Further information on the 
selection criteria are given in Sect. O The effect of choosing 
these lifetimes on the determination of the velocity distribu- 
tion functions /(|W|) in the magnitude bins are discussed 
in Sect, g^l 




[Fe/H]= -0.8, -0.68, -0.5, -0.28, 0., 0.20 
M/M,,,„=0.8: model A 
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Figure 2. The lower panel shows the age dependence of the 
absolute luminosity My for stars with M = 0.8 Mq (averaged 
over 0.1 Mq). Thin lines are the evolution for different metallici- 
ties [Fe/H] = -0.8; -0.68; -0.5; -0.28; 0.0; 0.20 and the thick line 
shows the present day luminosity of the stars as a function of age 
taking into account the age-metallicity relation. The upper panels 
show the luminosity of mass bins along the main sequence as a 
function of age taking the chemical enrichment into account. Ver- 
tical lines give the adopted lifetimes for each luminosity interval. 



3 OBSERVATIONAL DATA 

In the model we assume that each stellar sub-population is 
in dynamical equilibrium in the vertical direction. For the 
high luminosity end of the main sequence with average age 
below 1 Gyr this is not guaranteed. On the other hand the 
observed velocity distribution functions and also the verti- 
cal de nsity profile of A stars measured bv lHolmberg fc Flvnnl 
l|2000t) are hints for equilibrium distributions. Additionally 
we use large radii of 200 pc for these sub-samples and deter- 
mine a spatially averaged distribution. All sub-samples can 
be used as independent samples to determine the total ver- 
tical gravitational potential (/'(z). The velocity distribution 
functions of the samples depend on the age distri- 

bution of the stars. Therefore we use main sequence stars, 
where is a function of the lifetime of the stars. In or- 

der to determine the vertical velocity distribution function 
/(|H^|) in the solar neighbourhood we need kinematically 
unbiased samples with space velocities. 



For the determination of the AVR we use only the mea- 
sured velocity dispersions in the magnitude bins along the 
main sequence using stellar lifetimes. An indepen dent stellar 
sub-s ample are the McCormick K and M dwarfs l|Vvssotskvl 
Detected by a spectroscopic survey they are free from 
kinematic bias. Altogether 516 stars show reliable distances - 
almost all from the Hipparcos Catalogue - and space velocity 
components. All stars within the 25 pc sphere are also used 
to determine independently the velocity distribution func- 
tion /(I W^l) of stars with lifetime larger t han 12 Gyr (see Fig. 
[6|. For a sub-sample of about 300 stars IWilson fc WooUevI 
(|1970l ) estimated the Call emission intensity at the H and K 
lines in a relative scale allowing to construct six different age 
bins under the assumption of a constant SFR (see Jahreif5 
and Wielen, 1997). For each bin the vertical velocity disper- 
sion is determined. We do not use these sub-samples in our 
model for two reasons. The determination of the mean ages 
of the bins are inconsistent with our disc model and for the 
velocity distribution functions in the age bins the number 
of stars is too low to obtain reliable f{\W\){t). The AVR of 
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Figure 3. The HR-diagram shows all Hipparcos stars (grey dots) 
with Ut^ j-K ^0.15 within the distance limits chosen (see column 3 
in Table 1) The selected main sequence stars are over-plotted by 
larger dark dots. The eight magnitude bins are indicated by the 
horizontal dotted lines and the inclined line marks the selection 
including the turnofT stars. The full lines show the ZAMS with 
present day metallicity and the 'tracks' from Fig. [2] used for the 
selection process. 

this sample is shown in the lower panel of Fig. |8] only for 
comparison. 

For the determination of /(|W|) along the main se- 
quence we use the Hipparcos stars with good space velocities 
supplemented at the faint end down to My = 9.5 mag by 
stars from an updated version of the Catalogue of Nearby 
stars (CNS4) jjahreifi fc Wielen|[l997l '). i.e. also most of the 
CNS4 data rely on Hipparcos results. For the determina- 
tion of the space velocities good distances, proper motions 
and radial velocities are required. This was achieved in com- 
bining the Hipparcos data with radial velocities originating 
from an unpublished compilation of the "best" radial ve- 
locities for the nearby stars that was then extended to all 
Hipparcos stars. 

We apply a 2-step selection process. In the first step the 
absolute visual magnitude My of the Hipparcos stars was 
determined from the visual magnitude and the trigonometric 
parallax given in the Hipparcos catalogue. Only stars with 
a,r/7r smaller than 15% were taken into account. In the case 
of binaries resolved by Hipparcos the My of the brighter 
component was used. It was calculated from the combined 
magnitude and the magnitude difference measured by Hip- 



Table 1. Complete samples of nearby main sequence stars: Col- 
umn 1 lists the source catalogue. Column 2, 3 and 4 list the range 
in absolute visual magnitude, the selected distance limit, and the 
total number of stars available, respectively. Column 5 and 6 list 
the number of stars removed due to poor parallaxes (o-^/vr > .15) 
and unknown or poor radial velocities. In column 7 the remaining 
number of main sequence stars with good space velocity compo- 
nents is given and the last column contains the corresponding 
number of stars including the turnofT stars. 



source My cJhtti N n* no RV ^fin Nto 
[mag] [pc] 



Hip 


-2, -1, 


200 


304 


39 


15 


250 


931 


Hip 


1 


100 


242 





1 


241 


497 


Hip 


2 


75 


401 


8 


7 


386 


650 


Hip 


3 


50 


352 





4 


348 


491 


Hip 


4 


30 


172 


2 





170 


194 


Hip 


5 


25 


172 








172 




Hip 


6 


25 


170 





2 


168 




CNS 


7, 8, 9 


25 


525 




84 


441 





parcos. The stellar sample is divided into magnitude bins 
My = -1 ± 1.5, 1 ± 0.5, 6 ± 0.5, 8 ± 1.5 mag. In order to 
avoid a kinematic bias we restrict the distances of the stars 
in each magnitude bin to be well within the completeness 
limit of the Hipparcos catalogue determined by the magni- 
tude limit V ~ 7.3 mag of the Hipparcos Survey. The faintest 
magnitude bin is restricted to a distance of 25 pc relying on 
the completeness of the CNS4 (grey dots in Fig. |3]). 

In the second step we select in the HR-diagram 
(Fig. |3} a regime along the main sequence in or- 
der to exclude most of the turnoff stars as well as 
all giants and white dwarfs. As zero age main se- 
quence we used the {B — y. My ) -relation for the present- 
day metall icity z=0.025 originating from the Padua- 
isochrones ()Bertem et all Il994i : ICirardi et al.ll2002l 'l (URL: 
http://stev.oapd.inaf.it/~lgirardi/cgi-bin/cmd). Then all 
stars in the magnitude range My,ms±0.8 mag were selected. 
The reason for this additional selection is to match the mass 
tracks from stellar population synthesis up to the turnoff 
point described in Sect. 12.51 A selection of these tracks are 
over-plotted in the HR-diagram. For the brightest magni- 
tude bin the ZAMS is almost vertical and the stars redder 
than {B — V)zams + 0.10 were excluded. 

This selection is chosen to include the metal-poor MS 
stars below the ZAMS and to allow for a metallicity scatter 
of about 0.2dex as well as to include unresolved MS binaries 
above the single star MS. On the other hand most turnoff 
stars are excluded, because they would lead to an additional 
spread of stellar masses and lifetimes in the magnitude bins. 
For a consistency check we analysed a second set of stellar 
samples by including the turnoff stars bluer than the inclined 
line in (Fig. [3]) and using the corresponding larger lifetimes 
derived from stellar evolution. The properties of the result- 
ing samples of main sequence stars are collected in Table 

m 

For the determination of the velocity distribution func- 
tions /(jWI) we correct for the peculiar motion of the Sun 
us ing the local st andard of rest value Wq = 7km/s given 
m iDelhavd (|l965h that is in close agreem ent to the modern 
Wq_ = 7.17 ± 0.38 km/s determined by iDehnen fc Binneyl 
d 19981 ) from the Hipparcos data. The resulting normalised 
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distribution functions are shown in Fig. |6] with a binning of 
5km/s in \ W\. For the determination of the velocity disper- 
sions we exclude stars with high velocities (\W\ > 40km/s 
for Mv < 3.5 and \ W\ > 60 fcm/s for My > 3.5). 



4 PROPERTIES OF THE DISC 

We discuss first the algorithm to determine the parameters 
of best fitting models. Then we describe our fiducial model 
A in great detail and compare the properties with a selection 
of other models of similar fitting quality. 



4.1 Best fit procedure 

The disc model depends on a set of input functions and pa- 
rameters with strong correlations and degeneracies. Thus a 
simultaneous best fit of all parameters and an automatic se- 
lection of the astrophysical best model is impossible. There- 
fore it is very useful to sort the input quantities hierarchi- 
cally by there impact on the velocity distribution functions 
of MS stars and the structure of the disc. 

The fitting procedure can be split into three parts: 1) 
pre-processing to calculate a normalised model; 2) best fit 
of and scaling of the model; 3) post-processing to 

calculate the output and predictions (see figure |4]). We start 
with a reasonable set of mass fractions Qs.g.t of the com- 
ponents including the kinematics of the gas, thick disc and 
DM halo and fix the IMF. From large sets of template func- 
tions for the normalised SFR and AVR we select up to ten 
of each. We select for each SFR a metal enrichment law. 
Then we construct models for all (SFR, AVR) combinations 
by calculating the mass loss due to stellar evolution, solving 
the Poisson equation and determining the local stellar age 
distribution. With pre-calculated MS lifetimes we determine 
the normalised velocity distribution functions /i(|M^j) for all 
nine stellar sub-samples. 

In each model we determine the common scaling fac- 
tor (Je for all /i(|V K|) by a nonlinear Levenberg-Marquardt 
best fit algorithm (jPress et al.lll99l ). It turned out that a 
standard fit is very sensitive to the choice of the velocity 
range in |iy| due to the sparse population in the high er ve- 
locity bins. Therefore we tested the statistics from ILucvI 
Hooo) which is a modified statistics leading to a standard 
distribution also for sparse data sets. This algorithm is very 
stable and we fix the velocity range to |Vy | < 40 km/s for the 
bright bins with My < 3.5 and 60 km/s for the faint bins. 
For 90 degrees of freedom a reduced ~ 1-2 corresponds 
to a 10% level of significance that the data are consistent 
with the model. An inspection of the individual reduced 
for the /i(|H^|) shows that the brightest bin with x^ > 3.0 is 
inconsistent with the model on the 99.5% level for all mod- 
els. This result confirms that these very young stars are not 
in dynamical equilibrium. Nevertheless we include this bin 
in our calculation in order not to ignore the low velocity 
dispersion in that bin. Tests by excluding the first bin and 
by fitting also the normalisations rii of the /i(|VK|) does not 
change the results significantly. 

For further investigations we consider only models with 
X^ < 1.2 which already leads to strong constraints in the al- 
lowed heating laws AVR. All further iterations are done for 
individual pairs (SFR,AVR). In the next step the calculated 



'8,g,t 



, gas, DM 



SFR, AVR 



IMF 



[Fe/H](t) 





population 
synthesis 


mass loss, lifetimes 




r Poisson eq. 



scale heights, densities 
age distributions 



local [Fe/H] 



gas infall, yield 



PDMF 



f(W), 



Or 



G dwarfs 



star numbers 



best fit 



0„ 



scaling 



SFR, AVR selection 



local densities 
surface densities 
scale heights 
kinematics 



properties of components 



output 



densities, ages 
kinematics 
metallicities 
star counts 



Figure 4. Flowchart of fitting procedure. For a description see 
text. 



G dwarf metallicity is compared to the observations and a 
new chemical enrichment law is constructed to give a reason- 
able fit. Corrections via stellar mass loss and MS lifetimes 
are small and the new model can be used to derive the con- 
version factor from the IMF and the local present day mass 
function (PDMF). Now the IMF can be determined from 
the observed PDMF which is done in the current phase only 
for comparison and not iterated. 

For the next iteration step we need to scale the model 
to physical quantities. The first of the two free scaling fac- 
tors is CTc which is already fixed by the best fit proce- 
dure. For the second parameter we decided to fix the lo- 
cal stellar density including the thick disc contribution to 
PbO + pto = 0.039 Mq pc~^, which is the best observed quan- 



Milky Way disc model: Part 1 9 



tity for the local disc model (jjahreifi fc Wielerj[l997l ). We 
can rewrite equation |4] to derive Ze by 

2c = , — with Coff = /loff / 2c (25) 

47rGCoff Pso 

Then h^s can be scaled and the total surface density is de- 
termined by 

= (26) 

All other quantities can be scaled in a similar way. In this 
final iteration the input parameters for the surface densi- 
ties and the kinematics are adapted in such a way as to 
reproduce the known constraints on surface densities, local 
densities and scale heights. The thick disc parameters can 
be freely chosen in a wide range with very little influence on 
the other disc properties. The kinematic parameters of gas 
and DM halo are forced to get the correct gas scale height 
Zg ~ 100 pc and DM velocity dispersion an ~ 140km/s. 
The value of Qh and as a consequence the local halo density 
Pho are implicitly determined by fixing the other Q values. 
The value of Qg is chosen in a simple step to be consistent 
with the observed surface density of Eg « 10 — 13Mqpc~^ 
ijDickev fc Lockmanlll990h . Adapting the thin disc param- 
eter Qs is more sophisticated, because it has a significant 
influence on the shape of the gravitational potential. For 
the boundary conditions we have a relatively large freedom 
and compare models which fulfil different constraints on the 
thin disc scale height z^ or on surface densities. 

We discuss more details of parameter choice in the next 
subsections and a more detailed investigation of the best fit 
procedure and model selection will be given in Paper II. 

In the next subsections we discuss the selection of the 
main input functions SFR and AVR and other properties of 
the disc models including predictions for future observations. 

4.2 Star formation history and dynamical heating 

The main input functions that are to be determined, are the 
SFR and AVR. Since the SFR is a priori not strongly con- 
strained, we start with template sets of smooth functions of 
different analytic type: pure exponentials with different disc 
age and decay timescale; linear increase combined with ex- 
ponential decline; algebraic with linear increase and polyno- 
mial decline. For the heating function AVR we use templates 
of power laws with different power law indices and offsets in 
time to allow a nonzero present day value. A second set has 
a constant velo city dispersion fo r ages larger than 6 or 9 Gyr 
as proposed bv iFreemanI (jl99lh . A third set of polynomials 
was chosen to produce an extra flattening or steepening in 
the AVR around intermediate age of ~ 6 Gyr. 

In a preliminary study we select those parameter 
regimes which result in a reasonable < 1-2 for at least 
one heating function AVR. We found that the polynomial 
functions do not improve the fits for any SFR significantly 
compared to corresponding power law AVRs. The range of 
consistent SFRs, based on the local kinematic data alone, is 
weakly constraint: We find for the present day star forma- 
tion rate SFRp > 1.2MQ/pc^/Gyr (smoothed over the last 
500 Myr) and for stars older than 10 Gyr a contribution less 
than 25%. This means that the overall decline of the SFR is 
limited to a factor of 4-5. 



In the main step we restrict the models further to 
(SFR, AVR) pairs with < l-l- All models with x^ be- 
low this limit are consistent with the /i(|VK|) of the local 
MS stars. For each selected SFR the possible AVRs are re- 
stricted to at most a few. Since there is still quite a variety of 
very different models left over, we restricted the further in- 
vestigation to a representative selection of (SFR,AVR) com- 
binations. The fine-tuning of the model parameters (metal 
enrichment, gas and DM properties, surface density con- 
straints) is done only for these models. It turned out that 
these parameter adjustments do not change the best fit sig- 
nificantly (similar x^ and cto). In the models we fix the sur- 
face density of the gas to the observed value and the stellar 
surface density is determined by the local stellar density 
and the best fit CTc. Therefore the DM surface density, i.e. 
the mass fraction Qh, is mainly determined by the surface 
density of all components S (|z| < l.lkpc) ~ 75M|^ /pc^ 
iKuiiken fc Gilmord Il99ll ). (|Holmberg fc Flvnnl |2004| ). It 
must be determined by an iteration process, because chang- 
ing the Q's also changes the scale heights of all compo- 
nents. The second boundary condition is the scale height 
Zb of the thin disc. It is connected to Qa via the scaling 
length Zc (equ ation 1251). The standard v alue of 325±50pc 
determined by iBahcall fc Soneiral l| 19841 ) from star counts 
is still reliable and widely used. Since this scale height is 
calculated from a one component exponential thin disc, it 
corresponds to the mean exponential falloff in the range of 
\z\ ~ 100 — 1000 pc which is much closer to the mid-plane 
than our value Zs measuring the exponential scale height in 
the range \z\ ~ 600 — 1500 pc. An inspection of the density 
profile shapes (figure [TT|) shows that closer to the mid-plane 
the local scale height is ~20% larger than z^. Therefore we 
try to reach in the models z^ ~ 270 — 300 pc which is not 
possible in all cases. Reducing the DM fraction increases the 
thin disc fraction Qs and as a consequence the scale height 
Zb- As a bottom line there is a trade-off to reach both con- 
straints given by E(|2;| < 1.1 kpc) and z^. A general overview 
over all models and a detailed discu ssion of the paramete r 
dependence will be given in Paper II (|just fc Jahreifil|2009l ). 

We present four models A-D with very different SFR. 
For clarity we use model A as the fiducial model that will 
be described in great detail. The SFR of this model is ver y 
close to that investigated also bv I Aumer fc Binnev I l|2009l ). 
We discuss the similarities and differences of models B-D 
with respect to model A. The data of models E-G are added 
to show the effect of some parameter variations. The basic 
properties of the models are 

Model A: Fiducial model with maximum SFR 10 Gyr ago 
and mass fractions chosen to get simultaneously a small thin 
disc scale height of ~ 270 pc and a surface density E(j2j < 
l.lkpc) ~ 75Mq/pc^ 

Model B: Model with smallest reduced x^ = 0.96 of all 
models; optimised as model A 

Model C: Model with constant SFR; DM fraction in- 
creased to the upper limit in order to get a reasonable thin 
disc scale height of z^ ~ 300 pc 

Model D: Model with reduced disc age of t-p = 10 Gyr and 
exponential SFR with shortest consistent decline timescale. 

Model E: Model with same SFR, AVR as model A but 
neglecting the thick disc component. 
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For the dynamical heating lunction AVR we use a power law 

(30) 



6 4 
age [Gyr] 

Figure 5. The upper panel shows the SFRs of models A-D. The 
horizontal lines are the corresponding mean values (SFR) . The 
lower panel gives the AVRs of models A— D(x'^ values in paren- 
theses). Age is running backwards. 



Model F: Like model A but with larger local stellar den- 
sity. 

Model G: Like model F but with larger gas fraction in or- 
der to increase the local density further. 

The SFRs and AVRs of models A-D are shown in figure [S] 
together with the mean SFR determined by (SFR) = So/tp. 
A hst of model parameters and derived physical quantities 
discussed in this paper is given in Table [5] together with the 
corresponding data of other authors. 

For models A-D we give explicitly the input functions 
to allow the reader to make comparison calculations. We 
find for models A and B 



SFR(t) = (SFR) 

A : to 

B : to 

for model C 

SFR(t) = (SFR) 

and for model D 

SFR(t) = (SFR) 

D : = 8 



{t + to)tl 



with 



{t^ + tlY^ 

5.6 Gyr ; ti = 8.2 Gyr ; t„ = 9.9 Gyr 
1.13 Gyr ; t^ = 7.8 Gyr ; t„ = 11.7 Gyr 



(27) 



tp exp(-t/to) 
tc I — exp{—tp/tc) 
Gyr ; = 10 Gyr . 



(28) 



(29) 







( "'" ''"'^ 1 with 
\tp + roJ 


A : 


a 


= 0.375 ; TO = 0.17 Gyr ; tp = 12 Gyr 


B : 


Q 


= 0.5 ; TO = 0.5 Gyr ; tp = 12 Gyr 


C : 


a 


= 0.5 ; TO = 0.5 Gyr ; tp = 12 Gyr 


D : 


a 


= 0.5 ; To = 0.32 Gyr ; tp = 10 Gyr . 



Model C with constant SFR is of special interest, be- 
cause it is still used in many applications directly or implic- 
itly. Model C fits well the local kinematic data. Due to the 
large fraction of young and dynamically cool stars a larger 
initial velocity dispersion and a higher fraction of DM mat- 
ter is needed. The relatively large values for E(|z| < 1.1 kpc) 
and Za are still within the limits. 



4.3 Kinematics of the disc 

All presented models show reasonable fits to the set of nor- 
malised velocity distribution functions /i(|VF|) for the nine 
samples of MS stars described in section O because the to- 
tal reduced < l-l- For models A-D the total values 
are given in the lower panel of figure [5] In figure |6] the data 
are compared to model A in all samples. For comparison in 
each column one of the other models B-D are shown. We 
also computed the Xi of the individual distribution func- 
tions /idWI) to check the individual contributions. The val- 
ues are given in parentheses in figure [6] The distribution of 
the xf are consistent with a random distribution with one 
exception. In all models the brightest bin has an individual 
reduced Xi > 3 for eight degrees of freedom (= number of 
bins). The inset in the upper left panel of figure |S] shows the 
the distribution function of models A-D in log-scale for the 
bins between 10 and 20 km/s, where the significant devia- 
tions occur. The 1 — a standard Poisson noise of the data 
(shown by the errorbars) are not a precise measure of the 
contribution to x^ , since we use the statistics (for details 
see Paper II). This is a hint to a significant departure from 
the assumed equilibrium distribution. This is not surprising, 
since the lifetime of these stars is of the order of the vertical 
oscillation frequency in the disc and dynamical equilibrium 
cannot be achieved in this short time. In Fig. [7] we show 
for model A that the distribution functions /(|W|) includ- 
ing the turnoff stars agree just as well with the theoretical 
distribution functions using the appropriate lifetimes. 

The AVR can be observed only by kinematically un- 
biased stellar sub-samples with direct age determinations. 
In the lower panel of figure |8] the AVRs of the models A- 
D are reproduced from figure [5] in log-scale and compared 
to two data sets. The circles are the McCormick K and M 
dwarfs with ages determined from H and K line strength 
(Jahrcifi & Wielcn 1997) and the asterisks represent the F 
and G stars of GCS2 with good age determinations. The ages 
of the McCormick stars were derived by adopting a constant 
SFR and should therefore be consistent with model C. Cor- 
rections for a declining SFR as in the other models would 
shift the mean ages of the bins systematically to higher val- 
ues. It is interesting to note that the best AVR for models 
B and C are the same. The best fit scaling a a is different. 
This shows that the variation of the SFR has only a small 
influence on the shape of the AVR. The velocity dispersions 
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Table 2. Disc properties for models A— G and some models from the literature. Global parameters: value of the fit; maximum height 
2max, local densitypoi surface densities below 0.35 and l.lkpc, total surface density of the disc. Thin disc properties: maximum velocity 
dispersion a^, mean and present day SFR, local density pso; surface density up to ^maxi (half-)thickness h^ff, exponential scale height 
2a; Thick disc properties: velocity dispersion crt, local density pto; surface density St up to Zmax, exponential scale height zt; power law 
index at of sech™* fit. Gas component: local density pgo; surface density Eg up to Zmax, (half-)thicknoss ?tcfi,g; exponential scale heights 
Zg; DM halo: velocity dispersion (Tj,, local density pi^g; surface density Sjj up to Zmax- 



quantity 


unit 


A 


B 


C 


D 


E 


F 


G 


other sources 






fiducial 


min.x^ 


const. SFR 


lt.=10Gyr 


no. thick d. 


larger psO 


more gas 




red. X 




1.057 


0.960 


1.025 


0.994 


1.046 


1.056 


1.054 






kpc 


2.3 


2.7 


2.7 


3.1 


2.3 


2.3 


2.3 




PO 


M0/pc3 


0.088 


0.091 


0.092 


0.092 


0.088 


0.094 


0.101 


0.076\ 0.0982 


S(< 0.35) 


M0/pc2 


41.1 


41.3 


42.3 


40.7 


41.0 


43.5 


44.5 


4l3 


S« 1.1) 




73.5 


73.2 


79.3 


71.8 


72.5 


74.9 


74.5 


74±63, 71±6-* 


^disc 




45.2 


45.7 


42.1 


46.7 


41.8 


49.2 


50.4 


56±63, 48±9* 


(SFR) 


M0/pc2/Gyr 


3.75 


3.60 


3.25 


4.38 


3.99 


4.09 


4.02 




SFRp 


" 


1.43 


1.80 


3.33 


2.20 


1.55 


1.58 


1.42 




O-o 


km/s 


25.0 


28.4 


29.9 


28.0 


25.1 


25.0 


25.1 


17.5^ 


PsO 


M0/pc3 


0.037 


0.037 


0.037 


0.037 


0.039 


0.041 


0.041 


0.045\ 0.0442 


Ss 


M0/pc2 


29.4 


28.6 


26.3 


29.2 


31.4 


32.1 


31.6 


34.42 


hcB 


pc 


400 


389 


357 


398 


402 


386 


379 




2s 




274 


295 


303 


309 


277 


266 


265 




o-t 


km/s 


45.1 


45.4 


44.9 


44.9 




45.1 


45.1 


37^ 
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Figure 7. The normalised velocity distribution functions /(|H^|) 
of the four brightest bins as in Fig. |6] but with turnoff stars in- 
cluded. 



of the sub-samples of MS stars with no age subdivision are 
determined by the weighted mean over the lifetime with the 
local age distribution. 

For the AVRs of the different models we find power laws 



with indices 0.375 < a < 0.5. Th is is in the range of the 
classi cal value of 0.5 l)Wielenlll977^■ of .53 (Holmbc rg et al.l 
l2009h . and of 0.45 derived bv lAumer fc^Binnev I l|2009l) . The 
best fit slope depends strongly on the zero point for newly 
born stars. A higher initial velocity dispersion aw,o leads to 
a shallower heating function described by a larger power law 
index. All our models require a large (Jw,o and a steep rise 
of the AVR to reproduce the distribution functions of stars 
with lifetimes 0.5-3 Gyr. The maximum velocity dispersion 
of the oldest thin disc stars is 25-30 km/s dependent on the 
SFR. 

The upper panel of Fig.|8]shows the excellent agreement 
of the model with the data from the nearby stars. The black 
triangles are the velocity dispersions in the magnitude bins 
of the main sequence stars with lifetimes from Fig. (2] The 
full circle is the velocity dispersion of the McCormick stars. 
The underlying distribution functions of these data sets were 
used for the best fitting. The grey triangles are the velocity 
dispersions including the turnoff stars in the five magnitude 
bins with Mv < 4.5 mag. As expected, they show a system- 
atically larger velocity dispersion and have larger lifetimes 
compared to the pure IVIS samples. These data confirm the 
underlying assumption of a continuous disc heating and dy- 
namical equilibrium of the stellar sub-samples. 
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Figure 6. The normalised velocity distribution functions for the different sub-samples given in Tablc[T]and for the McCormick 

stars (lower right). Symbols are the data and full lines represent the fiducial model A. In each column the distribution functions of one of 
the comparison models are added: B: best x^; C: const. SFR; D: exp.SFR and lifetime=10 Gyr. The inset in the upper left panel shows 
the part of (in log-scale), where the significant deviations occur. 



Since the shape of /(|W|) depends on the age distri- 
bution of the stars, it depends also on the height z above 
the mid-plane. Fig. |9] shows the variation above the plane 
for stars with lifetimes larger than 12 Gyr (upper panel) and 
stars with a lifetime of 6 Gyr (lower panel; corresponding to 
magnitudes A/v = 3.5 — 4 mag). The vertical gradient in the 
velocity dispersion of all thin disc stars is slightly stronger 
than that of the stars with lifetimes larger than 12 Gyr, be- 
cause the dip near the mid-plane is enhanced due to the 
contribution of stars with smaller lifetimes. The profile is 



shown in figure [TO] for models A and C. The large fraction 
of young stars in model C with constant SFR leads to a 
significantly larger gradient. 



4.4 Density profiles 

The shape of the vertical density profiles of all components 
and sub-populations are determined by their kinematics and 
the self-consistent gravitational potential. A general feature 
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Figure 8. Upper panel: The symbols show the velocity disper- 
sion of the stars along the main sequence (magnitude bins) as a 
function of mean lifetime determined from stellar evolution (Fig. 
[2]|. The full line gives the result of the fiducial model A. The other 
lines are for the comparison models B— D. Lower panel: The lines 
show the velocity dispersion as a function of age (AVR) for the 
stellar sub-populations of the models. Full symbols a re the age 
group s of the McCormick K and M dwarfs (iJahreifi &: WielenI 
Il997t) and the asterisks are from GCS2. 



of the vertical density profiles of the stars and of the gas is 
the flattening at the galactic plane. All profiles are between 
an exponential profile and that of an isolated isothermal pro- 
file given by a sech^ function. In the upper panel of figure fTUl 
the profiles of model A are shown. The stellar disc profile of 
model C (also plotted in figure [T0)| is significantly flatter and 
deviates from model A in the regime \z\ ^ 500 pc and above 
\z\ ~ 1500 pc. The second plot in figure [TOl shows the differ- 
ence in the force that is a measure of the total surface 
density up to \z\, of models A and C. Below 500 pc the pro- 
files are very similar despite the differences in the mass frac- 
tions an d the kinematics. Th e corre spon ding surface density 
values o f lKuiiken fc Gilmord (Il99ll ) and lHolmberg fc Flvnnl 
(|2004l ) are added. 

For the gas, the thickness is ~150pc compared to 
the exponential scale height of ~100pc. The density pro- 
file and the surf ace density are consiste nt with the ob- 
served Hl-profile (|Dickev fc Lockman|[l990l ) : 0.014 M0/pc^ 
5.0Mq/pc^, 177 pc for the mid-plane density, surface den- 
sity, thickness) adding about 50% of H2 with smaller scale 




Figure 9. The normalised velocity distribution functions /(|W|) 
(for dw = 5km/s) at different heights z for model A. The upper 
panel shows of stars with lifetimes larger than 12 Gyr including 
the thick disc contribution, whereas the lower panel is for thin 
disc stars with lifetime 6.0 Gyr. 



height l|Bronfman et al.lll988l ) and applying the correction 
factor of 1.4 for Hehum and heavy elements. 

Due to the gravitational potential of the disc, the local 
density of the dark matter halo is ~50% larger than the DM 
density at Zmax and ~20% larger than the mean DM density 
given by Eh/(2z 

max I • 

A measure of the profile fiattening is given by the ratio 
of the thickness and the exponential scale height. The half- 
thickness of MS stars hms as function of lifetime is shown in 
the lower panel of Fig.[T2]for all models A-D. The differences 
above a lifetime of 6 Gyr is significant leading to differences 
in the determination of surface densities from local densities 
as used in the estimation of the SFR and IMF from local 
data. For model A additionally the exponential scale height 
Zoxp as function of lifetime and the half-thickness /iago as 
function of age are plotted. For young populations Zcxp is 
much smaller than hms, because the density profiles have 
Gaussian and not exponential wings. For longer lifetimes 
/ims is ~50% larger than Zcxp- 

The density profiles of MS stars differ in shape from 
the profiles of the sub-populations of single ages. The lower 
panel of Fig. [TT] shows the normalised density profiles with 
lifetimes according to the samples used for the model. They 
are steeper than the corresponding profiles of the sub- 
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Figure 10. Vertical profiles of the model A and of the thin disc 
of model C (SFR=const.). From top to bottom: Density profiles; 

force complemented by the corresponding surface den sities of 
iKuiiken fc Gilmore) l ll99l[) and lHolmbere fc Flvnnl | |2004) : veloc- 
ity dispersion of the thin stellar disc; V-band luminosity weighted 
mean metallicity of the thin stellar disc. 



populations with the same age and significantly shallower 
than the density profile using the mean age of the sub- 
population. This difference is quantified by the thicknesses 
ftage and ftms showu in the lower panel of Fig. 1121 

An independent test of the model is the comparison 
with dir ectly observed dens it y pro files of MS stars as de- 
rived by iHolmberg fc Flvnnl ( 2000l ) for two sub-samples of 
A stars (with < My < 1 mag) and early F stars (with 
1 < Mv < 2.5 mag). The mean stellar lifetime in these mag- 
nitude bins are 0.3 Gyr and 1.0 Gyr, respectively, including 
turnoff stars. In order to minimise the systematic asymme- 
try of the profiles (mainly for the A stars) we applied an 
additional offset of the solar position in z of -1-5 pc. The re- 
sult is shown in figure 1131 For the A stars all profiles are 
shifted by a factor of two to avoid overlaps. Over-plotted 
are the normalised profiles of models A-D and F,G. The ob- 
served profile for the F stars is in very good agreement with 
all models. The observed profile of the A stars is slightly 
steeper than the model and requires a larger local density 
or an even shorter lifetime. Models F and G show the effect 
of a larger local density in two steps (see table [2}. Since we 
do not expect that these stars are already in dynamical equi- 
librium, the deviations are in a reasonable range. In the Be- 
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Figure 11. Upper panel: Normalised density profiles of the stellar 
sub-populations as a function of age. Lower panel: Same for the 
main sequence stars as a function of lifetime. The set of lifetimes 
corresponds to the mean lifetimes of the sub-samples used for the 
velocity distribution functions (see Fig. |6]l. 



sangon model (|Robin et al.ll2003 ) Einasto laws are used for 
the density distribution of isothermal sub-populations. The 
thicknesses (=fiattening in their model) are determined by 
solving the Poisson equation. The shapes of the density pro- 
files deviate systematically from our results and this means 
that dynamical equilibrium is achieved only approximately. 
A detailed investigation of the differences to our model is 
beyond the scope of this paper. 



4.5 Age distributions 

The age distributions of the stellar sub-samples of MS stars 
depend on the lifetime and on the vertical structure. The 
upper panel of figure [141 shows for model A the normalised 
local age distribution for different lifetimes. The compari- 
son with the SFR shows clearly the over-representation of 
young stars in the solar neighbourhood due to the increasing 
thickness with age. The vertical dilution is quantified by the 
(half-)thickness /ims (lower panel of Fig. I12p . 

A comparison of the local age distributions of stars with 
lifetime larger than the age of the disc for models A-D are 
shown in the lower panel of figure [T^ In model A the local 
age distribution of stars with lifetime larger than the disc age 
varies by less than a factor of two around the mean value. 
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Figure 12. Top panel: Comparison of the density profiles for thin 
disc stars with lifetimes larger than the age of the disc. Bottom 
panel: The (half-)thickness of the sub-populations as a function 
of age /lage and of the MS stars as a function of lifetime /ims- For 
comparison the exponential scale heights Zexp for the MS stars 
are included. The horizontal lines are the overall thickness /left 
and the exponential scale height Zs of all stars. 



Figure 14. The upper panel shows the normalised local age dis- 
tributions of MS stars with different lifetimes. The middle panel 
shows the age distribution of stars with lifetime larger than 12 Gyr 
at different heights above the Galactic plane. The lower panel 
shows the local age distributions of models A-D. 
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Figure 13. D ensity profiles for the A and F star sample from 
iHolmberg fc Fl vnn (2000) with an arbitrary vertical offset (full 
and open symbols below and above the mid-plane, respectively. 
The lines are models A— D and F,G for corresponding mean life- 
times 0.4 Gyr and 1.0 Gyr for the A and F stars, respectively. The 
profiles for the A stars are shifted additionally by a factor of two. 



iBinnev et aD (|200Ci ) and ICignoni et all ()2006l ) propose a 
constant local age distribution in this sense. This is con- 
sistent with our model taking the uncertainties of isochrone 
ages into account. In model B there is a lack of old stars 
which cannot be tested by direct measurements, because for 
this age range the age determinations are very uncertain. 
Model C with constant SFR predicts a stro ng dominance 
of you ng s tars in contras t to t he findings of iBinnev et al.l 
(|200d ) and lCignoni et all l|2006l ). 

The age distributions are a strong function of z above 
the plane. The middle panel of Fig. [14] shows the lack of 
young stars in steps of Az = 300 pc above the mid-plane. 



4.6 Metallicity 

Since MS luminosities and lifetimes depend on metallicity, 
we include a simple analytic metal enrichment law [Fe/H](t). 
We adopt a generalised form of a closed box model with 
a Schmid t star formation law for the oxygen abundance 
[0/H](t) (|just et al.lll99^ ;) and transform it by a linear ap- 
proximation to [Fe/H] from lReddv et aP (|l993l ) for thin disc 
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stars, I.e. 

[Fe/H](t) 
[0/Hl(t) 

^oxy {t) 



2.67[0/H](t) 

lg(Zoxy(t)) 

+ (-^ox; 



(31) 

lg[l + ln(l + g(t/tp)'-] 
ln(l + q) 



There are four fitting parameters q, r, [Fe/H]Q, [Fe/H]p wliicli 
determine tlie scaling and the curvature. We tested a large 
range of parameters: initial metallicity [Fe/H]Q=-0.6. . .-0.8, 
present day metallicity [Fe/H]p=-0.1. . . +0.2, linear increase 
up to steep early enrichment covering [Fe/H]=-0.4. . . +0.1 
at an age of 6 Gyr. In order to reproduce the high metallic- 
ity tail without intrinsic scatter we investigated also some 
AMRs with an additional accelerated enrichment in the leist 
2 Gyr. In every case it turned out that the local G dwarf 
metallicity distribution is a strong restriction for the AMR 
for any given SFR and AVR. The parameters for models 



A-D 


are 










A : 


q = 2 


; r = 0.55; 


[Fe/H]„ = 


-0.6; [Fe/H]p = 


0.02 


B : 


q = 2 


; r = 0.55; 


[Fe/H]„ = 


-0.8; [Fe/H]p = 


0.01 


C : 


9 = 3 


; r = 0.70; 


[Fe/H]„ = 


-0.67; [Fe/H]p = 


= -0.01 


D : 


9 = 2 


; r = 0.55; 


[Fe/H]„ = 


-0.6; [Fe/H]p = 


0.01 



The metal enrichment laws are shown in the lower panel 
of figure rroi For model A [O/H] is also plotted. The local 
metallicity distribution is very sensitive to the present day 
metallicity due to the shallow slope. 

From the metal enrichment the local metallicity dis- 
tribution for stars with lifetime larger than the age of the 
disc is calculated. Before binning the theoretical distribu- 
tion we add a Gaussian scatter which represents an intrinsic 
scatter and observational errors. We tested a range up to 
0.2 dex and found a best value of 0.13 dex corresponding 
to a FWIIM=0.31 dex in order to reproduce both the nar- 
row width of the maximum and the high metallicity wing 
in the observations (see figure I15|l . The local metallicity 
distribution is determined from the Copenhagen F and G 
star sample (|Nordstrom et al. .2004l l selecting all stars with 
masses 0.84 ^ M/Mq ^ 0.90 up to the completeness limit 
r < 40 pc. The lower mass limit of 0.84Mq was chosen in 
order to avoid incompleteness for the less luminous metal 
rich stars. The upper limit was set to be consistent with a 
lifetime larger than 12 Gyr for the selected stars. 

For model A the comparison of the derived local metal- 
licity distribution for stars with adopted lifetime larger than 
12 Gyr and the data of the selected mass bin is shown in the 
upper panel of Fig. 1151 The results for models B-D are very 
similar. The middle panel of Fig. [15] shows the predicted 
metallicity distribution in model A for stars with lifetime 
6 Gyr compared to the observed distribution with lifetime 
larger than 12 Gyr in the solar neighbourhood. The lower 
panel shows the same comparison but at \z\ = 900 pc. 

The metal enrichment [O/H] and the SFR can be em- 
bedded in a local chemical evolution model with primordial 
gas infall. Since the gas infall rate is an additional free func- 
tion, the metal enrichment cannot be derived directly from 
the SFR, IMF and stellar yields. We proceed in the follow- 
ing way. For the oxygen enrichment we adopt instantaneous 
recycling and mixing to determine the infall rate of primor- 
dial gas. The mass loss from stellar evolution is taken into 
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Figure 15. Upper panel: The metallicity distribution of model A 
compared to the observations from GCSl in the stellar mass range 
0.84 ^ M/Mq ^ 0.90. Lower panels: The metallicity distribution 
for dwarfs with lifetime 6.0 Gyr in comparison to the dwarfs with 
lifetime > 12 Gyr in the solar neighbourhood and 900 pc above 
the mid-plane. 



account. Then the oxygen yield in solar units is given by 



(Zoxy) + Zo:^y,p-^ + Zoxy,0-^ (32) 
JO JO 

Here (Zoxy) is the mean metallicity of all born stars (not 
in log-scale and weighted by mass). Go is the initial sur- 
face density of gas. Eg the present day value and So the 
total amount of born stars. We start with a negligible ini- 
tial amount of gas Go and find for the yields in solar units 
of models A-D j/oxy = 1.1135, 1.1301, 1.1290, 1.1350, respec- 
tively. The infall rate, the surface density of gas and stars 
and the metal enrichment laws are shown in Fig. 1161 

4.7 Thick disc 

The thick disc properties in the solar neighbourhood are not 
well determined. In many investigations a local stellar con- 
tribution of ~ 5% and a velocity dispersion of ~45 km/s are 
adopted. These values lead to a dominance of thick disc stars 
compared to the thin disc at d istances z > 1 kpc which is 
directly observed for K dwarfs iPhleps et al.ll2005l ). But in 
a recent paper local contribu tions up to 20% an d metallici- 
ties up to solar are claimed (|Bensbv et al.ll2007l ). Since the 
method of the latter paper relies on a statistical separation 
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Figure 16. Upper panel: Gas infall rates for models A— D. Middle 
panel: Cumulative stellar and gas surface density (thick and thin 
lines, respectively). Lower panel: Metal enrichment laws [Fe/H] 
of models A— D. The horizontal thin black line at solar metal- 
licity is for guiding the eye. Additionally the oxygen enrichment 
[O/H]=0.375[Fe/H] is shown for model A (pink dotted Une). 



from the kinematics by adopting Gaussian velocity distri- 
bution functions for the thin and for the thick disc, it un- 
derestimates strongly the thin disc contribution at the high 
velocity tail. In our models we use standard values for the 
thick disc and investigate the influence of the thick disc on 
the other model properties. We split the local stellar density 
into a thin disc and a thick disc component. For models A-D 
we set at ~ 45km/s and choose Et leading to a local stellar 
mass contribution of ~ 5%. The local mass fraction of the 
thick disc equals the thin disc density at z = 1.1 — 1.25 kpc 
in the models. We adopt an age of 12 Gyr for the thick disc 
leading to two consequences. Firstly the high mass end of 
the MS is missing and the fraction in number of low mass 
stars is a factor of 1.84 larger than the mass fraction. This 
has to be taken into account when comparing to results of 
star counts. For model A the number of K dwarfs in the thick 
disc equals that in the thin disc at z = 950 pc very c lose t o 
the observed K dwarf density profile of lPhleps et al.l (|2005h . 
Secondly for the kinematics /(|W^|) the thick disc contribu- 
tion is only taken into account for stars with lifetime larger 
than 12 Gyr. Here the correction factor 1.84 is also included. 

In all cases the density profile of the thick disc can be 
fitted by a sech" {z / aht) profile to better than 3%. The cor- 



responding best fit coefficients are given in tabled For the 
thin disc profile a corresponding fit deviates up to 25%. 

We have calculated two alternative thick discs in model 
A with smaller and larger velocity dispersion, respectively. 
The surface density of the thick discs are adjusted to yield 
the same crossing point of the thin and thick disc profile 
(upper panel of figure [17)) . the local density contributions 
of the thick discs are 5.7%, 8.6% and 4.9% for discs 1,2, 
and 3, respectively. The contribution of the thick disc to the 
velocity distribution function is shown in the lower panel 
of Fig. 1171 All these models are still consistent with the 
data in the -range of 40-60 km/s. From the distribution 
functions we can conclude that a kinematically distinct thick 
disc cannot exceed a local density contribution of 10% in 
mass density. 

Model E is a model without thick disc component but 
with the same AVR and SFR as in model A. For the min- 
imum no alteration in the AVR and SFR are needed. 
The best fit (Je is slightly higher. The total star formation 
5*0 and the densities of the thin disc are larger containing 
now all stars in the solar neighbourhood. The corrections to 
all other parameters are within a few percent. The effect of 
the thick disc on the velocity distribution function of stars 
with lifetime larger than 12 Gyr is shown in the lower panel 
of Fig. 1171 Only the high velocity tail is affected. 



4.8 Initial mass function 

We have fixed the IMF for the best fitting process (equation 
I23p . We can use the local luminosity function to test the 
adopted IMF. For the determination of the IMF from the lo- 
cal luminosity function, the luminosity function is first con- 
verted to a mass fun ction dN = f(M)dM using the trans- 
formation formulae o f lHenrv fc Mc Carthvl l| 19931 ). corrected 
m iHenrv et al.] (|l999t ). a nd extended to bright stars accord- 
ing to Schmi dt-Kaler in ISchmidt-Kalerl (| 19821 ) , which was 
confirmed by I Andersen I (|l99ll ). The star numbers are nor- 
malised to a sphere with radius R=20 pc. The result is shown 
in the l ower panel of Fig. 1181 and compared to the PDMF 
given in lKroupa eFahl (|l993 ) (KTG93). In our model we use 
the system mass function excluding resolved B components 
of binaries. In addition turnoff stars were excluded to be 
consistent with the main sequence lifetimes used in the dy- 
namical model. Therefore our PDMF is systematically be- 
low KTG93. The PDMF is also corrected for the vertical 
gradient of the density profiles in the observed volume. In 
the brightest bin with lifetime 0.15 Gyr the density in the 
mid-plane is 50% larger than the mean density in the 200 pc 
sphere. 

In order to determine the conversion factor from the 
local PDMF to the IMF due to the finite lifetime and the 
vertical thickness hms, the main sequence lifetime is needed. 
We use the lifetimes estimated from the stellar evolution 
tracks shown in Fig. [5] for the different magnitude bins, 
where the metal enrichment and the relative weighting due 
to the increasing number of stars with decreasing mass in the 
mass interval is taken into account. These data are shown 
in the upper right panel of Fig. [18] with a comparison of 
the lifetimes determ ined directly from evolutionary tracks of 
Girardi et al .l |2004l) and with the analytic fitting formula of 
Eggleton et al. (1989). Any systematic variation of the life- 



times result in significant changes of the conversion factor 
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Figure 17. Upper panel: Density profiles of models with differ- 
ent isothermal thick disc components: Thin disc stars of model E 
without thick disc and model A with different thick discs; thick 
disc parameters {rTt,ci,ht) are indicated. Lower panel: Reproduc- 
tion of f\W\ from the last two plots of figure [B] but in log-scale 
to demonstrate the effect of the thick disc contribution. Thick 
lines are the combined distribution functions and thin lines are 
the thick discs only. 



and therefore the IMF. The conversion factors are a combi- 
nation of the dilution by the increasing thickness with age 
(measured by /ims) and the fraction of born stars still on the 
main sequence due to the finite lifetime. The last factor de- 
pends strongly on the SFR. The upper left panel of Fig. 1181 
shows the conversion factor (=IMF/PDMF) as a function of 
lifetime at different heights above the mid-plane for model 
A. 

In the lower panel of Fig. [18] the PDMF and the IMFs 
of model A and C are shown with different symbols. The 
IMFs of model B and D are very similar to that of model 
A. We compared different best fits in log-scale to test the 
significance of different IMFs. The highest mass bin is not 
consistent with the Scalo IMF (Eq. [23}, but the dominating 
contribution to the rms value of 4.98 originates from the 
mass range 0.4 < M/M^T) < 1. A si milar best fit of a Kroupa 
IMF KTG93 ([Kroupa et al.lll993t ) leads to a rms value of 
5.67 for a = -2.7 at M/Mq > 1.0. If we allow the slope of 
the Scalo IMF to vary at the high mass end for M /Mq > 
2.0 we find a best fit value of a = —4.16 and a rms of 
2.81. A similar best fit of a Kroupa-like IMF as used in 
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Figure 18. Upper left panel: Conversion factor (=IMF/PDMF) 
for model A that must be applied to the PDMF to derive the 
normalised IMF from local star counts. Upper right panel: Stel- 
lar lifetimes as function of stellar mass to relate the conversion 
factors at the left to the PDMF below. Dots are the values used 
in the models (from Fi g. [2ll and lines are val ues determined by 
iGirardi et al] ||2004|) and lEegleton et al.l lll989l) . Lower panel: Full 
circles show the PDMF of our sample and crosses the PDMF of 
KTG93. The IMF from models A and C are given by the plus 
signs and the triangles, respectively. The best fit IMF and corre- 
sponding PDMF for model A (full red lines) and model C (dotted 
blue lines) are over-plotted. The dashed black lines show the Scalo 
IMF and IMF for comparison. 



lAumer fc Binnev I (|2009h leads to q = -3.06 at M/Mq > 
1.0 with rms=5.96. 

An inspection of the shape of the IMFs of models A 
and C shows that the break near 1 Mq seems artificial. The 
break is obvious in the PDMF, but the reason for this is that 
the correction factor due to the finite lifetime starts to apply 
above ~ IMq . Therefore we determined an alternative IMF 
by fitting power laws in two mass regimes only, where we 
fixed the break mass Mq by minimising the rms of the best 
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fit. The best fits are 



5 SUMMARY 



diV/MQ = 
A: a = 

C : a = 



Nq {M/MqY d (^/Mq) with 



(33) 



-1.46 ±0.10 
-4.16 ±0.12 

No = 272/Mq at 

-1.70 ±0.13 
-4.07 ±0.11 

TVo = 388/Mq at 



0.08 ^ M/Mq < 1.72 
1.72 s: M/Mq < 6.0 

Mo/Mq = 1.72 , 

0.08 ^ M/Mq < 1.35 
1.35 ^ M/Mq < 6.0 

Mo/Mq = 1.35, 



where A^o gives the normalisation in the 20 pc sphere. The 
slopes at the high mass end are very similar for models A-D, 
because the IMFs are essentially shifted vertically relative 
to each other. The best fits of the Kroupa-like IMFs show 
that forcing the slope in the mass range 1 — 1.5Mq to be the 
same as at the high mass end results in a flatter IMF and 
a much larger rms value. At that stage of the investigation 
we abstain from final conclusions, because we did not dis- 
cuss possible biases in the PDMF from local star counts. We 
expect that the feedback of a corrected IMF with a steeper 
slope at the high mass end to the disc model via the mass 
loss is very small. Additionally it affects only the velocity 
distribution functions for lifetimes smaller than 10* yr sig- 
nificantly. 

The strongest constraints on the IMF at high masses 
and the present day SFR is the observed number of A and 
late B stars in the solar neighbourhood. Since the bright 
stars with My < 0.5 mag are observed in a sphere with a 
radius of 200 pc, the sample size is a direct measure of the 
local surface density. Therefore the conversion from the IMF 
to the mean SFR in the last few 100 Myr depends only on 
the lifetime of the stars. That means, a higher present day 
SFR requires a shorter lifetime for the A stars or a steeper 
IMF. 

As an ex ample the constan t SFR used in the Be- 
sangon model (|Robin et all |2003| ) relies mainly on the as- 
sumption of a steep IM F in the stellar mass range 1-3 Mq 
l|Havwood et alj 1 1997^1. Thi s mod el and also the SFR de- 
termination of Vergelv et al.l (|2002l ') show that the slope of 
the IMF in the mass range of 1-3 Mq that covers roughly 
the lifetime range of 0.2-10 Gyr, is crucial for the deriva- 
tion of the long-term SFR from star counts. A steep IMF 
in that mass range introduces a strong bias to young ages, 
because the observed number of stars above 2 Mq predict 
too many young lower mass stars via the steep IMF. This ef- 
fect is strengthened by the fitting procedure, because direct 
age determinations by isochrone fitting are less significant 
in the lower mass range. The result is a bias to young ages 
in the stellar mass range covering the lifetime above 1 Gyr. 
An underestimation of the SFR for ages above 1 Gyr is the 
consequence. 

The positi on of the Sun is probably ab out 20 pc above 
the mid-plane (|Humphrevs fc Larsenlll995l ). For the deter- 
mination of the mid-plane density this offset can play a 
significant role for sub-populations with small scale height 
(My < 1.5 mag). We corrected for that implicitly, since we 
determined the mid-plane density of these magnitude bins 
for the PDMF using spheres with large radii (see Table [1]), 
where the offset can be neglected. 



We presented a new disc model for the thin disc in the solar 
cylinder based on a continuous star formation history (SFR) 
and a continuous dynamical heating law (AVR) of the stellar 
sub-populations. This new model combines and improves the 
advantages of several different attempts to model the verti- 
cal structure in the solar neighbourhood: We used a sequence 
of isothermal sub-popula tion s in dynamical equ i libriu m as 
iBahcall fc Soneiral (|l984') and lAumer fc Binnev I l|2009D did; 
We used the full velocity distribution functions f\W\ as 
iHolmberg fc Flvnnl l|2000l ) did, and not only the velocity 
dispersions (the AVR); We solved the Poisson equation self- 
consistently including the thick disc, gas and DM halo con - 
tribution as done in the Besangon model (|Robin et al.ll2003l ). 
A chemical evolution model with reasonable gas infall rate, 
which was tuned to reproduce the local [Fe/H] distribution 
of G dwarfs, is included. This enables us to apply correct 
MS luminosities and lifetimes. Additionally our model is in- 
sensitive to the IMF, because it is based on the normalised 
velocity distribution functions of MS stars. 

We determined pairs of (SFR, AVR) by a best fit of the 
local kinematics. The SFRs which are consistent with the 
MS velocity distribution functions show a decline factor be- 
low five down to unity (= const. SFR). The strongest feature 
that would distinguish between a constant and a declining 
SFR is a direct determination of the age distribution of low 
mass stars in the solar neighbourhood. 

Despite the large variety of SFRs there is a strong cor- 
relation to the AVR. For each SFR the slope and maxi- 
mum velocity dispersion of the AVR are well determined. 
For the AVR we find a power law with indices between 0.375 
and 0.5. The range of models is consistent with the results 



of Binnev et al.l (I2OO OI) for the local age distribution and 
lAumer fc Binnev I I2OO9) for the SFR. 

Applying the stellar lifetimes and the new scale height 
corrections to the PDMF results in an IMF that shows only 
one break point at 1.7Mq and a steep falloff at high masses. 

The density profile of an isothermal thick disc compo- 
nent can be fitted very good by a sech" profile, where at de- 
pends on the velocity dispersion. The most prominent effect 
of thick disc stars is the enhancement of the high velocity 
wings of K and M dwarfs in the range of 40-60 km/s. From 
that we can exclude a heavy thick disc with distinct kine- 
matics and more than 10% contribution to the local stellar 
density. Changing the thick disc parameters leads to slight 
variations of the thin disc properties (mainly by assigning 
part of the stellar disc to the thick disc), but has a negligible 
influence on the normalised SFR and AVR of the thin disc. 

A variety of predictions can be made from the new disc 
model. The density profiles of MS star sub-populations de- 
pend on the lifetime of the stars and are significantly differ- 
ent to density profiles of single age sub-populations. The 
shape is neither exponential nor sech'^ and can be char- 
acterised by the (half-)thickness and the exponential scale 
height. From the vertical density profiles MS stars number 
densities as function of colour and apparent magnitude can 
be predicted. Applying these number densities with observed 
'Hess' diagrams from large surveys like the catalogues of the 
Sloan Digital Sky Survey SEGUE/SDSS enables us to re- 
strict the parameters of the SFR further. We determined 
vertical gradients in the kinematics which will be tested with 



20 A. Just & H. Jahreifi 



Radial Velocity Experiment (RAVE) data. Age and metal- 
licity distributions of stellar sub-populations as a function 
of 2 above the galactic plane are predicted. 

The future plan is to extend the local model to a com- 
plete disc model of the Milky Way that provides a fully 
self-consistent connection of stellar densities and kinematics. 
Ultimately this kind of detailed model is essential to under- 
stand the large data sets as already available from SDSS and 
which are expected on a much higher level in amount and 
precision by PanSTARRS and the astrometric Gaia satellite 
mission. 
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